Discrete soliton collisions in a waveguide array with saturable nonlinearity 
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We study the symmetric collisions of two mobile breathers/solitons in a model for coupled wave 
guides with a saturable nonlinearity. The saturability allows the existence of breathers with high 
power. Three main regimes are observed: breather fusion, breather reflection and breather creation. 
The last regime seems to be exclusive of systems with a saturable nonlinearity, and has been previ- 
ously observed in continuous models. In some cases a "symmetry breaking" can be observed, which 
we show to be an numerical artifact. 
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I. INTRODUCTION 

Since the 1960's, a great number of papers have consid- 
ered the properties of solitons in nonhnear optic media 
with a Kerr-type (cubic) nonlinearity. This media can 
be modelled by the cubic Nonhnear Schrodinger (NLS) 
Equation. As it is well known, the NLS equation is in- 
tegrable and, in consequence, solitons interact elastically 
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More recently, several authors have studied the prop- 
erties of solitons in photo- refractive media 0. In this 
case, the equation describing these media is a modifica- 
tion of the original NLS, which consists in substituting 
the Kerr nonlinearity term by another one of saturable 
type. This Saturable (SNLS) Equation is nonintegrable 
and the soliton collision processes are inelastic, leading 
to annihilation, fusion or creation of solitons This 
last phenomena consists of the appearance of three soli- 
tons after the collision of only two of them. Another 
important feature of the SNLS is that the behaviour of 
the solutions is quite generic, being independent of the 
details of the mathematical model. 

The discrete version of the NLS equation can be used 
to describe nonlinear waveguide arrays within the tight 
binding approximation The existence and proper- 
ties of mobile discrete breathers/solitons in DNLS lat- 
tices has been considered in a number of studies (We use 
the terms breathers and solitons interchangeably in this 
context, also intrinsic localized modes). An early brief 
study showed that breathers could propagate along 
the lattice with a small loss of energy, and could be- 
come trapped by inhomogeneities in the lattice. Later, a 
more detailed study suggested that "exact" travelling 



breathers might exist, at least for some parameter ranges. 
The reviews 0, ^ refer to many other papers in this area. 
More recently, work has concentrated on breathers with 
infinite oscillating tails 0, although the question of the 
existence of exact breather solutions which tend to zero 
as n — *■ ±cxD has not yet been resolved. Given the long 
history of mobile breather solutions of this equation, it is 
rather surprising that a systematic study of the collision 
of two breathers in the DNLS model has only recently 
been carried out (We mention also that collisions 

have been studied in generalised nearly integrable DNLS 
model in [Hill). 

Recently, some studies have considered the existence of 
mobile breathers in waveguide arrays in photo-refractive 
crystal, described by a DNLS equation with saturable 
nonlinearity In particular, these papers consid- 

ered a discrete version of the Vinetskii-Kukhtarev model 
0, The key difference between the cubic DNLS 

equation and the saturable DNLS equation is that in the 
later, the Peierls-Nabarro barrier (the energy difference 
between a bond-centred and a site-centred breather with 
the same power) is bounded and, in most cases, smaller 
than in the former . It allows the existence of mobile 
breathers of high power. 

The aim of the present paper is to study breather- 
breather collisions in a saturable DNLS equation and to 
compare the results with those obtained in the continuous 
SNLS and the discrete cubic equation. 



II. NUMERICAL RESULTS 

The system we consider is governed by the following 
equation of motion 
u 

■ /3 , , 1" I, + iu„+l - 2Un + Un-l) = 0. (1) 
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This model has two conserved quantities: the Hamilto- 
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nian H = E„[/31og(l 
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p] and the 
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power (or norm) P = J2n I'^nP- 

In order to reduce the dimension of the large parameter 
space to be considered, we have fixed (3 to (3 — 2. Higher 
values of (3 lead to solutions that only can be moved for 
a restricted set of power values ld\. 

A moving breather Vn (t) is obtained by adding a thrust 
5 to a stationary breather u„, so that: 



Vn{0) = u„exp(i qn). 



(2) 



Notice that this procedure of obtaining moving breathers 
is similar to the marginal mode method introduced in 
for Klein-Gordon lattices. 

In the following, we consider the collision of two iden- 
tical breathers moving in opposite directions with the 
same thrust q. Analogously to Ref. 0|, we consider 
both inter-site (IS) and on-site (OS) collisions. 

The collision scenario we observe for small P is quite 
simple: there exists a critical value qc below which 
breathers form a bound state, and above which, breathers 
are reflected (See Fig. ^-b for examples of these two 
cases). It can be observed that the bound state "oscil- 
lates" after the collision. The amplitude of these oscil- 
lations decreases when approaching to the critical point, 
whereas their "period" increases. (Note that the "reflec- 
tion" case could equally be regarded as a transmission 
case as the two breathers are indistinguishable. In the 
case of reflection/transmission, there is some loss of en- 
ergy of the two breathers). 

For high values of P, the above scenario takes place, 
except that, for high values of q, breather creation is ob- 
served. Figure^ shows an example of such a collision. 
This behaviour is similar to the soliton creation observed 
in the saturable continuous models and will be analyzed 
in more detail below. The different regimes in the (P, q) 
plane are depicted in Figure |21 for both IS and OS colli- 
sions. Furthermore, Fig.^shows the values of the critical 
value of q separating merging and reflecting regimes, as 
a function of the power P. It can be seen that, for most 
choices of P, both values are close. This is different from 
the cubic DNLS case where the critical values of the 
OS case is an order of magnitude higher than the ones 
for the IS case. The likely explanation is that in the sat- 
urable case, the PN barrier is small (for our choice (3 = 2, 
the absolute value of the barrier is smaller than 0.01). 

For high values of P, we have also observed the merging 
of two breathers with symmetry breaking, as reported in 
[TOi] . This symmetry breaking manifests as a movement 
of the final bound state to left or to the right accom- 
panied by the appearance of a total lattice momentum, 
defined hy p = i Y,n ( 1 V'^ - 1 V'n ) • Since the equa- 
tion, the initial conditions, and the boundary conditions 
are symmetric, this state must be a numerical artifact, 
as suggested in To test this hypothesis further, we 
performed some runs with either (a) increased numerical 
accuracy in the numerical integration routines, or (b), the 
addition of some very small random noise to the initial 
conditions. In case (a), the onset of symmetry breaking 
is shifted to longer times, whereas in case (b), symmetry 




FIG. 1: Typical power density plots for (a) bound state for- 
mation (P = 10, q = 0.1), (b) reflection (P = 10, g = 0.2), 
and (c) breather creation (P = 70, q = 0.5). In all cases, OS 
collisions are considered, although these pictures do not vary 
considerably for IS collisions. 



breaking is observed at shorter times. These numerical 
results confirm that symmetry breaking is a numerical 
artifact caused by random rounding errors breaking the 
symmetry of the problem. However these "spurious" re- 
sults are interesting in their own right as they suggest 
that at these higher values of P, the stationary breather 
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formed after collision is more easily set into motion by 
a very small perturbation. To check that the other phe- 
nomena we observe is not due to numerical artifacts, we 
have carried out similar tests on other runs showing dif- 
ferent phenomena. No such sensitivity to random errors 
of the accuracy of the integrators is observed. 




FIG. 2: Different regimes observed in (a) IS and (b) OS colli- 
sions. The colours represent the following: white-merge to a 
single breather; black-reflection; and red-breather creation 
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FIG. 3: Critical value of the initial thrust q for (a) IS and (b) 
OS collisions versus P. 



phenomenon is similar to the soliton bistability observed 
for SNLS solitons in It consists of the existence 

of a minimum in the dependence of the soliton width 
with respect to the peak intensity. This dependence is 
monotonically decreasing in the cubic NLS, and thus the 
soliton in a saturable medium has a Kerr behaviour for 
small peak intensities (or power). In the discrete case, as 
the width is less well-defined, we have considered instead 
W = |uiP/|moP, where n = is chosen as the centre (or 
peak) of the breather. Fig.Elshows W versus P display- 
ing a similar behaviour as in the continuous case. 

The analysis given above also explain why the results 
of l^l (i.e. only merge and reflection regimes take place) 
are found for small values of the power. We note also 
that this creation process may be related to the phenom- 
ena of the fission of a coupled two-breather state into a 
stationary and a moving breather in the DNLS equation 

la. 



(a) 



(b) 
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FIG. 4: Minimum value of the power density at the collision 
point after collision as a function of P and q. Left (right) panel 
corresponds to inter-(on-) site collisions. We have supposed 
that the power trapped in the trapping regime is zero in order 
to clarify the figure. 
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III. BREATHER CREATION 

We proceed to analyze the breather creation process, 
as it is most noteworthy phenomenon that appears in 
the saturable case in comparison to the cubic one. From 
Figs. 13 and 01 we can conclude that the conditions for 
breather creation are that P and q are above a threshold 
value. 

This result can be explained with the aid of Fig. [Sj 
where the density power of the collision point, for the 
cases of reflection and breather creation, is shown. It 
can be seen that the power density oscillates after the 
collision, and its minimum is zero for the case of no cre- 
ation. The minimum power density after the collision 
for all the simulations (neglecting the trapping regime) 
is represented in Fig. ^ In consequence, the trapped 
power must be above a threshold so that breather cre- 
ation occurs. It can be explained by the fact that, for 
a stationary breather to have a "saturable" behaviour, 
its power should be higher than a threshold value. This 



FIG. 5; Time evolution of the power density at the collision 
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The left panel corresponds to a reflection case 



{q — 0.3, P — 70) and the right panel to a creation case 
(g = 0.7, P = 70). 



IV. EFFECT OF PHASE IN BREATHER 
COLLISIONS 

To complete the paper, we give a brief study of 
the effect of considering a phase difference between the 
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(a) (b) 
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FIG. 6: Representation of the "breather width" W (defined as 
W = |«ip/|itoP) versus the power of a stationary site-centred 
(left) and a bond-centred (right) breather. 



The extra power available to breathers in the SDNLS 
case results in the new phenomena of breather creation 
in a discrete model. Additionally, the scenario in the sat- 
urable DNLS case seems to be much "cleaner" than in the 
(a) (b) 
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breathers, in a similar fashion to Ref. This is 

achieved by introducing a factor exp(i^) in one of the 
breathers. Fig. |7| shows the final amplitudes Ai^2 and 
velocities V1.2 of both breathers as a function of (p. 




FIG. 7: Final velocities (a) and amplitudes (b) with respect 
to the phase for breathers with P = 20 and q — 0.25 and an 
OS collision. 



It appears that the final velocities are smooth functions 
of (j), showing a strong phase effect, with the V2 curve 
following the Vi curve, phase-shifted by n. The V values 
vary from around 0.2 to 0.8. The amplitude dependence, 
on the other hand, is much smaller but shows a much 
more irregular behaviour as a function of (j). Clearly the 
discreteness of the lattice is featuring strongly here in 
this latter case. 

Finally, in Fig. |S1 we show the relation between the 
outgoing velocities as varies through 2tt, analogously 
to Fig. 3 of Here the relatively smooth behaviour 
over a large range of V values is clearly shown. 



FIG. 8: Relation between the final velocities for breathers 
with q = 0.25 and (a) P = 20; (b) P = 10. OS collisions are 
considered in both cases 



cubic DNLS case on a coarse scale, with a strong but sim- 
pler threshold effect. These facts may be an advantage 
in some applications, such as multi-port optical switch- 
ing. There are still a number of details in the fine-scale 
structure which are as yet unexplained. These may per- 
haps be understood through the application of a future 
variational study. It would also be interesting to extend 
this study to consider the collision of two non-identical 
breathers. 
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Notes added in proof Since the original version of this 
paper was written, Khare et al. |23| have published an 
exact stationary breather solution for an equation closely 
related to namely 



1 + A^IV'ril 



rV'n + [i^n+l - '^i'n + ^n^l) = 0. (3) 



Solutions of (j3Jl can be mapped into solutions of iQJ by 
the (invertible) transformation 

■ipn{t) = — exp{wt//i} Un{t), (3 = vl[l, 



V. CONCLUSIONS 

We have analyzed the coUisional behaviour in a sat- 
urable DNLS model, finding close analogies to the con- 
tinuous NLS equation. Breathers can merge, refiect or be 
created (although breather annihilation is not observed). 



so the solutions given |20| can be mapped into solutions 
of lU . However the localized stationary breather solution 
of ,20] only exists for /3 > 2, and hence are not relevant 
to our discussions which focus on the /3 = 2 case. It 
would be interesting to extending the calculations in this 
paper to other values of /3 to see if the presence of these 
stationary solutions affected the results given here. 
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